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1. INTRODUCTION 

The fundamental inverse problem for vowel articulation is concerned [1-5] with the 
determination of the geometry of the human vocal tract from some data. In this paper, 
we consider various types of scattering data in the frequency domain resulting from a 
unit amplitude, monochromatic, sinusoidal volume velocity sent from the glottis, and we 
analyze whether each data set uniquely determines the shape of the vocal tract, or else, 
what additional information may be used for the unique recovery. 

Let us use x to denote the distance from the glottis and / for the length of the vocal 
tract. Hence, the lips are located at a distance / from the glottis. Typically, / varies 
between 14 cm and 20 cm, usually smaller for children than for adults and smaller for 
females than for males [1,2,5]. Even though the vocal tract is not a right cylinder, to a 
good approximation it can be treated as one [3,4]. 

We will let A(x) denote the cross sectional area as a function of the distance from 
the glottis, and we suppose that A(x) is positive on [0, /]. Assuming that the propagation 
is lossless and planar (these assumptions are known [3,4] to be reasonable), the acoustics 
in the vocal tract is governed [1-5] by the first-order linear system of partial differential 
equations 



where t is the time variable, the subscripts x and t denote the respective partial derivatives, 
\i is the air density, c is the speed of sound, v(x,t) is the volume velocity of the air flow, 
and p(x,t) is the pressure at location x and at time t. 

The volume velocity is equal to the product of the cross sectional area with the 
average velocity of the air molecules crossing that area. The pressure is the force per 
unit cross sectional area and is exerted by the moving air molecules. The air density at 
room temperature is fx = 1.2 x 10~ 3 gm/cm 3 . The speed of sound varies slightly with 
temperature, and c = 3.43 x 10 4 cm/sec in air at room temperature. In our analysis of 




(1.1) 
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the inverse problem, we assume that the values of fx and c are already known. There is no 
loss of generality to start the time at t = 0. 

By using v x t = Vt x , we can eliminate v in (1.1) and obtain Webster's horn equation 
' [A(x)p x (x,t)] x - \ Ptt (x,t) = 0, xe(0,l), t>0. 



Letting 

$(x,t):=y/A(x)p(x,t), (1.2) 
we find that $ satisfies the plasma-wave equation 

®xx(x,t) - \$tt(x t t) =Q(x)$(x,t), xe(0,l), t>0, (1.3) 

where we have defined 



with the prime denoting the x-derivative. The quantity Q is called the relative concavity 
of the vocal tract or the potential. Separating the variables as 

&(x,t) := ij(k,x)e lkct , (1.5) 

we find that ip(k, x) satisfies the Schrodinger equation 

ip"(k,x) + k 2 il>(k,x) = Q(x)il>(k,x), xe(0,l). (1.6) 

kc 

The frequency z/ is measured in Hertz and related to k as v = — . Informally, we can 
refer to k as the frequency even though the proper term for k is the angular wavenumber. 

In order to recover A, we will consider various types of data for k G R + resulting 
from the glottal volume velocity v(0,t) given in (4.1). As our data sets, we consider the 
absolute value of the impedance at the lips, the absolute value of the pressure measured at 
a microphone placed at some distance from the lips, the absolute value of the pressure at 
the lips, the absolute value of the transfer function from the glottis to the lips, the absolute 

3 



value of the impedance at the glottis, the absolute value of a Green's function for (1.3) 
measured at the lips, and the real or imaginary part of the reflectance at the glottis. 

The inverse problem of recovery of A can be analyzed either as an inverse spectral 
problem or as an inverse scattering problem. In the inverse-spectral formulation, in addi- 
tion to a boundary condition at the the glottis such as (2.1), a boundary condition is also 
imposed at the lips. The imposition of the boundary conditions at both ends of the vocal 
tract results in standing waves that are related to an infinite sequence of discrete frequen- 
cies. It was established by Borg [6] that Q can be recovered by using two such infinite 
sequences of discrete frequencies corresponding to two sets of boundary conditions. It then 
follows [3,4,7-11] that A can be recovered from two infinite sequences of constants. For 
example, such sequences can be chosen as the zeros and poles [7,8] of the input impedance 
or the poles and residues [9] of the input impedance. 

In the the inverse-scattering formulation, a boundary condition is imposed at only one 
end of the vocal tract — either at the glottis or at the lips. Then, the measurement of the 
acoustic data used in the recovery of A is performed at the same end or at the opposite 
end. If the boundary condition and the measurement occur at the same end of the tract, 
the corresponding inverse problem is usually known as a reflection problem. On the other 
hand, if the boundary condition and the measurement occur at different ends, then we 
have a transmission problem. The methods based on the inverse scattering formulation 
may be applied either in the time domain or in the frequency domain, where the data set is 
a function of t in the former case and of k in the latter. We refer the reader to [3,4,12-15] 
for some approaches as time-domain reflection problems and to [16] for an approach as 
a time-domain transmission problem. Our approach in this paper is a frequency-domain 
approach, where the analysis in Sections 6-10 may be viewed as that for a transmission 
problem and the analysis in Sections 11 and 12 may be viewed as that for a reflection 
problem. 

Our paper is organized as follows. In Section 2 we review some preliminary material 
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related to the Schrodinger equation and introduce the selfadjoint boundary condition in- 
volving cot a given in (2.1), the Jost solution /, the Jost function F a , and the scattering 
coefficients T, L, and R. In Section 3 we briefly review the recovery of Q, cot a, F a , /, T, L, 
and R from the data {\F a (k)\ : he R + }. In Section 4 we obtain some explicit expressions 
for the pressure and the volume velocity in the vocal tract in terms of A, /, and F a , and 
we also show that cot a appearing in (2.1) is directly related to the physical parameters 
A(0) and A'(0). In Section 5 we introduce the relative area [r](x)] 2 and express it in terms 
of the Jost solution, cot a, and the scattering coefficients. In Sections 6-12 we analyze the 
recovery of Q, and A from various data sets. The data set used in Section 6 is the 
absolute value of the output impedance at the lips. In Section 7 it includes the absolute 
value of the pressure measured at a microphone placed at some distance from the lips. In 
Section 8 it is the absolute value of the pressure measured at the lips. In Section 9 the 
data set includes the absolute value of the transfer function from the glottis to the lips. In 
Section 10 it is the absolute value of an analog of the Green's function introduced in [17] 
for (1.3), in Section 11 the absolute value of the input impedance at the glottis, and in 
Section 12 the real or imaginary part of the reflectance. Finally, in Section 13 we present 
some examples to illustrate the theoretical results presented in the earlier sections. 

2. PRELIMINARIES 

In this section we review the scattering data related to the potential Q appearing 
in the Schrodinger equation on the half line R + with the selfadjoint boundary condition 
[18-21] 

sina- (p'(k,0) + cos a • </?(&;, 0) =0, (2.1) 

where a is a number in the interval (0, it) identifying the boundary condition at x = 0. We 
can relate the half-line Schrodinger equation to (1.6) by assuming that Q(x) = for x > I. 
Note that the mapping a i— > cot a is one-to-one and from (0, n) onto R. 

Let / denote the Jost solution [18-23] to the half-line Schrodinger equation. It is 
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uniquely determined by the asymptotic conditions 

f(k,x) = e ikx [l + o(l)], f'(k,x) = ike ikx [l + o(l)], x -> +00. 
Since Q vanishes when x > /, we have 

f(k,l) = e ikl , f'(k,l) = ike ikl . (2.2) 

The Jost function F a associated with the half-line Schrodinger equation with the 
boundary condition (2.1) is defined as [18-21] 

F a (k) := -i[f'(k, 0) + cot a ■ f(k, 0)]. (2.3) 

Let us emphasize that the subscript in F a identifies the boundary condition at x = and 
it does not indicate any partial derivative. It is known [18-21] that 

F a (-k) = -F a (k)*, keK, (2.4) 

where the asterisk denotes complex conjugation. 

We assume that Q is real valued and integrable on (0,1) and that there are no bound 
states for the half-line Schrodinger equation with the boundary condition (2.1). The ab- 
sence of bound states is equivalent [18-21] to assuming that F a (k) has no zeros on I+, 
where I + := z(0, +00) is the positive imaginary axis in the complex plane. It is known 
[19-21] that either F a (0) ^ or F a (k) has a simple zero at k = 0; the former is known as 
the generic case and the latter as the exceptional case. The exceptional case corresponds 
to the threshold where the number of bound states can be changed by one under a small 
perturbation of the potential. 

By using the extension Q(x) = when x < 0, we can relate f(k, 0) and f'(k, 0) to the 
scattering coefficients in the full-line Schrodinger equation. We have [19,20,22,24] 

/(M) = ^, /'(M) = ^, (2.5) 
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where T and L denote the transmission coefficient and the left reflection coefficient, re- 
spectively, associated with Q. The right reflection coefficient R is given by 

It is known [19,20,22,24] that 

T(-k)=T(k)*, R(-k) = R(k)*, L(-k) = L(k)*, keR. (2.7) 

The absence of bound states for the full-line Schrodinger equation is equivalent [19,20,22,24] 
for T(k) not to have any poles on I + , and this is also equivalent [25] for F^^ik) not to 
have any zeros on I + . 

3. RECOVERY OF Q FROM |F Q | 

In the absence of bound states, the fundamental inverse scattering problem for the 
half-line Schrodinger equation with the selfadjoint boundary condition (2.1) consists of 
determining Q and cot a from various types of scattering data. In this section we review 
the solution to this inverse problem when the data set is {|F a (/c)| : k G R+}. 

Theorem 3.1 Assume that Q is real valued, measurable, and integrable for x G (0,/). 
Then, the data set {|F a (/c)| : k G R + } uniquely determines Q(x) for x G (0,/) and cot a. 
The same data set also uniquely determines the corresponding Jost solution f(k,x) and 
the scattering coefficients T(/c), R(k), and L(k). 

Below we outline some steps involved in the solution to the inverse problem stated 
in Theorem 3.1. As seen from (2.4), \F a (k)\ is an even function of k G R, and hence the 
data sets {|F Q (/c)| : k G R + } and {|F Q (/c)| : k G R} are equivalent. By using the data 
{|F a (/c)| : k G R} as input in the Gel'fand-Levitan method [18-21], we form the kernel 
function 

i r°° r p 1 

— 1 (cos kx) (cos ky) , 



1 f°° 
G a {x,y) := - / dk 

7T ./_™ 



\F a (k)\* 
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and then solve the Gel'fand-Levitan integral equation 

h a (x,y) + G a (x,y) + / dz G a (y, z) h a (x, z) = 0, < y < x. (3.1) 

Jo 

The solution to (3.1) is known [19-21] to exist and to be unique. Once h a (x, y) is obtained, 
we recover the potential as 

d 

Q(x) = 2 — h a (x, x ), x E (0,1), 
ax 

where x~ indicates that the limit from the left must be used in the evaluation. We also 
recover the boundary condition as 

cot a = —h a (0, 0). 



Alternatively, we can proceed [21] as follows. Let 

Aa(fc):=-1 + ^^, fceC+ 
F a {k) 



(3.2) 



where we use C + for the upper half complex plane and C+ for C + U R. Then 

k 2 



Re[A a (fc)] = -l + 



k e R. 



(3.3) 



\F a (k)\ 2 ' 

From the data {\F a (k)\ : k £ R} we first construct the function A a (A;) via the Schwarz 
integral formula as 



K(k) 



dt 



r>i .1 -oc t — k — z0+ 



-1 + 



k e C+, 



(3.4) 



\F a (t)\ 2 _ 

where the quantity i0 + indicates that the values for real k should be obtained as limits 
from C + . Next, F a (k) is obtained from \F a (k)\ by using 

\og\t/F a (t)[ 



F a (k) = k exp 



V 7H J-c 



dt 



t - k - zo+ r 



k e C+. 



Then, we have 



1 



/(fc,0) = -F a (AO[l + A a (fc)], keC+, 



(3.5) 



(3.6) 



f'(k,0)=iF a (k) 



+ lim [kA a (k)} 

cota = — % lim [/cA a (/c)], 

k — >oo 



fc e C+. 



(3.7) 
(3.8) 



where the limit in (3.8) can be evaluated in any manner in C+. Having both f(k,0) and 
f'(k,0) in hand, we can construct all the quantities that are relevant in the scattering 
theory for the Schrodinger equation. For example, the scattering coefficients for the full- 
line Schrodinger equation can be obtained as 



T(k) 



2ik 



L(k) = 7 



ikf(k,0)-f(k,0) 



ikf(k,o) + f(k,oy v ' ikf(k,o) + f(k,oy 

R{k)= -i*/(-M)-/'(-M) 



(3.9) 



(3.10) 



ifc/(fc,0) + /'(M) ' 

Having obtained such quantities, we can construct the potential by using any one of the 
various methods available [19,20,22,24]. For example, we can use the Faddeev-Marchenko 
method [19,20,22,24] and get 

Q(x) = -2-^K(x,x + ), x E R, 
ax 

where K(x,y) is obtained by solving the (left) Faddeev-Marchenko integral equation 

/>oo 

K(x,y)+R(x + y)+ dz R(y + z) K(x, z) = 0, -oo < x < y, (3.11) 

J X 

with the kernel 

1 C°° 

R (y)--=2^J dkR(k)e lk y. 
The Jost solution f(k,x) can directly be obtained from K{x,y) as 



f(k,x) 



+ 



f 

J x 



dyK(x,y) e 



iky 



(3.12) 



Let us remark that that, in order to obtain {A a (k) : k £ R} from {|F Q (/c)| : k E R}, 
instead of using (3.4) we can equivalently construct the real and imaginary parts of A a (k) 
via (3.3) and 



Im[A„(fc)] 



zicpv^ dt 



7V 



t - k 



-1 + 



keK, 



where CPV indicates that the integral must be evaluated as a Cauchy principal value. 
Consequently, cot a can be recovered from the equivalent form of (3.8) given by 

cotct= lim (klm[A a (k)]) . 

fc^oo 

4. PRESSURE AND VOLUME VELOCITY 

When the vocal tract area function A is known, via (1.4) we can evaluate the potential 
<5, solve the corresponding Schrodinger equation, and obtain the Jost solution f(k,x). In 
this section, with the help f(k, x), we express the pressure and volume velocity correspond- 
ing to the glottal volume velocity 

v(0,t) = e lkc \ t>0. (4.1) 
It is known [19,20,22] that f{k, •) and f(—k, •) are linearly independent for each k G 



C+ \ {0}. Hence, the general solution to (1.6) can be written as a linear combination of 
f(k, •) and f(—k, •). From (1.2), (1.5), and (1.6), we see that the pressure has the form 

p(x,t) = P(k,x)e lkct , (4.2) 

with 

P(k, x) = —±= [a(k) f(-k, x) + b(k) f(k, x)] , (4.3) 

y/A{x) 

where a(k) and b(k) are coefficients to be determined. When x > /, the pressure p(x, t) 
should be a wave traveling outward from the lips and should not contain the part propor- 
tional to e tk ( ct+x ^ traveling into the mouth. Thus, with the help of (2.2) we see that we 
must have b(k) = in (4.3). Hence, (4.3) is reduced to 

P(k, x) = — L= a(k) f(-k, x), xe [0, /]. (4.4) 

y/A(x) 

Our aim is to determine the value of a(k) in terms of the pressure P(k, I + r) measured 
by a microphone placed at a radial distance r from the lips. The relationship between 
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the pressure measured at the microphone and the volume velocity at the lips is explicitly 
known and is given by (cf. (3.1) of [5]) 



p(l + r, t) 



4rrr 



v t (l,t- r/c), 



(4.5) 



where we recall that fi is the air density and c is the sound speed. 

From (4.1), (4.2), and the first line of (1.1), for the x-derivative of the pressure at the 



glottis we get 



P'(k,0) = 



ikcfi 

W) 



Note that from (4.4) through differentiation we obtain 

1 



P'(k,0) = 



a(k) 



/'(-M)-#^/(-M) 



2A(0) 



where we have used 



[VMx)Y A'(x) 



^Mxj 2A(x) 
A comparison of (4.7) with (2.3) shows that, by choosing 



cot a = 



A'(0) [y/A(x)]'\ x=0 



2A(0) 



we can write (4.7) as 



P'(k,0) 



a(k)F a (-k). 



Comparing (4.6) and (4.10) we get 

a(k) 



ckn 



^W)F a {-ky 

and hence we can write (4.4) in the equivalent form 

= cfikf{-k,x) 

Using (4.2) and (4.11) in the first line of (1.1), we get 



x e [o,t\. 



v t (x,t) = 



ckA(x)e %kct d 
y/A(Q)F a (-k) Tx 



f{-k,x) 
^A{x) 
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(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



xe[0,l], t>0. (4.12) 



ik + 



In particular, from (4.12) we obtain 

vAl, t — r/c) = . 

y/W) F a (-k) 

Using (2.2) in (4.13) we have 
Finally, comparing (4.5) and (4.14) we get 

ckfiy/A(t) e ikc(t-r/c-l/c) 

or equivalently, with the help of (4.2), we have 



t > r/c. (4.13) 



A' (I) 
2 A{1) 



t > r/c. 



ik + 



A' (I) 
2 A{1) 



t > r/c, 



F a (-k) 



4nr y/A(0)P(k,l + r) 



ik + 



A'(l) 
2A(l) 



From (4.16) we can conclude that k appears as ik in P(k, I + r) and hence 



(4.14) 



(4.15) 



(4.16) 



P(-k,l + r) = P(k,l + r)*, keR. 



(4.17) 



We emphasize that P(k, x) given in (4.11) is valid only when x G [0, /], and hence P(k, l + r) 
is not obtained from (4.11) by replacing x by / + r there. Finally, we remark that, with 
the help of (4.1) and (4.12), we obtain 



v(x, t) 



i ^fA(x)e ikct 
^W)F a {-k) 



/'(-M-^/C-M) 



x E [0,1], t > 0. (4.18) 



5. AREA AND RELATIVE AREA 

Let us view (1.4) as the zero-energy Schrodinger equation, and consider the initial- 
value problem 

y" = Q(y)y, x e (0,1), 

(5.1) 

y(0) = ^MOj, y'(0) = -y/MG) cot a, 
12 



where cot a is the quantity in (4.9). Let yi(x) and yz(x) be any two linearly independent 
solutions to (5.1) on the interval [0,/]. Then, the unique solution to (5.1) can be written 

as 



yi{x) y 2 (x) 
-1 yi (0) 2,2(0) 
cot a yl(0) y' 2 (0) 
where [F;G] := FG' — F'G denotes the Wronskian. Let us define 



y{x) = 

[yi(x);y 2 (x) 



(5.2) 



rj(x) 



(5.3) 



We will refer to [r](x)] 2 as the relative area of the vocal tract. Let us remark that we can 
write (5.3) in the equivalent form 



A(x) = A(0) [rj(x)f 



x e [0,/]. 



(5.4) 



Recall that the Wronskian of any two solutions to the Schrodinger equation is indepen- 
dent of x, and [yi(x); y2(x)] ^ if and only if y\ and yi are linearly independent on [0, /]. 
For example, we can choose y\ and y 2 as the zero-energy Jost solutions g\ (0, x) and g T (0, x), 
respectively, for the full-line Schrodinger equation where the potential agrees with Q(x) 
on the interval (0, /), is zero when x < 0, but is a real- valued, measurable, integrable func- 
tion with a finite first moment when x > I. Let r(/c), £(k), and p(k) be the corresponding 
transmission coefficient, the left reflection coefficient, and the right reflection coefficient, 
respectively. Generically, we have r(0) = or equivalently [t7i(0, x); g v {Q, x)] ^ 0. In the 
exceptional case, we have r(0) ^ or equivalently [(?i(0, x); g r (0, x)] = 0. 

In the generic case, using [19,20,22,24] 

2ik 



[g r (k,x);gi(k,x)] 



g r (k,x) 



r(kY 

gi(-k,x) +p(k)gi(k,x) 



(5.5) 



<7r(0,0) 



r(k) 
^(0,0) = 0, 
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we can write (5.2) as 



7]{x) 





2 

1 0i (0,0) 

-cot a (71(0,0) 

where the overdot denotes the /c-derivative 



-- f(0) yi(0, x) i <7i(0, x) - ^ p(0) yi(0, x) 



1 




(5.6) 



In the exceptional case, we can choose 

Vi{x) = gi(0,x), y 2 {x) = gi(0,x) nJn 

In this case, we have 

1 + ^(0) 



dz 

gi(o, z y 



yi(o) = 



-, y'i(0) = 0, y 2 (0) = 0, y' 2 (0) = 



r(0) 



r(0) ' ^ ^ yi (0) 1 + £(0Y 

with yi(0) ^ because [19,20,22,24] we have -1 < £(0) < 1. Hence, from (5.2) we get 

nx dz 



7](x) = gi(0, x) 



1 / 

Jo 

1 01(0,0) 

cot a 



gi(o,zy 
o 
l 



<7i(0,0) 



(5.7) 



Theorem 5.1 The relative area [r](x)] 2 for x G [0,/] is uniquely determined by the data 
{\F a (k)\ : k G R + }. Equivalently, n(x) for x G [0,/] is uniquely determined from {Q(x) : 
x G (0, 1), cot a}, where cot a is the constant in (4-9). 

PROOF: From Theorem 3.1 we know that {\F a (k)\ : k G R + } uniquely determines the 
potential Q(x) for x G (0,/) and the constant cot a. From (4.8), (4.9), and (5.3) we see 
that 77(0) = 1 and r/(0) = — cot a. Thus, rj(x) is uniquely obtained by solving the initial 
value problem in (5.1) in the special case A(0) = 1. | 

It is possible to construct rj(x) from {|F a (fc)| : k G R + } as follows. With the help 
of (3.4)-(3.12), we can construct the corresponding right reflection coefficient R{k), the 
transmission coefficient T(fc), and the Jost solution f(k, x). Hence, in the generic case every 
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term appearing on the right hand side of (5.6) can be constructed from {\F a (k) | : k G R + }, 
and we get 



r](x) 



-lT(0)/(0,x) if(0,x)- l -R(0)f(0,x) 



1 /(0,0) 
-cot a f(0,0) 

In the exceptional case, from (5.7) we get 



V(x) = f(0,x) 







1 



dz 



(5.8) 



1 /(0,0) 
1 



cot a 



/(0,0) 



f(k,x) = 



(5.9) 



Note that we can write the Jost solution f(k,x) as a linear combination of g\(k,x) 
and gi(—k, x), where g\{k, x) is the quantity appearing in (5.5). With the help of (2.2) we 
obtain 

gi(k,x) gi(-k,x) 

1 9\(kJ) g\(-kj) 
ik g[(kj) g[(-kj) 

We can also express F a (k) with the help of g\(k,x). To do so, we can obtain /(/c,0) and 
f'(k,0) from (5.9) and use (2.3) to get F a (k). Alternatively, by using [cf. (2.5)] 

01 (M) = Tm - > ^i(/c,0)=z/c 



ikl 



2ik 



r(k) 



r(k) 



we can write F a (k) with the help of the transmission and left reflection coefficients asso- 
ciated with g\(k, x). 



6. RECOVERY FROM THE IMPEDANCE AT THE LIPS 



The impedance at the lips is defined as 



Z(k,l) := 



p{l,t) 
v(l,t)' 



(6.1) 
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Using (4.2), (4.11), and (4.18) in (6.1), we get 



Z(k,l) 



2ick/i 



(6.2) 



2ikA(l)+A'(l)' 

Thus, we can only hope to get A(l) and A' (I) from Z(k, I). We can refer to Z(k, I) as the 
output impedance because the volume velocity is from the glottis to the lips. In Section 11 
we will analyze the impedance at the glottis, which we can identify as the input impedance. 

Note that from (6.2) we get 

4c 2 £;V 



|z(M)r = 



k G R. 



(6.3) 



4k 2 A(lf + A'{1) 2 ' 

By using (6.2) at two distinct real k values, say k\ and k 2 , we can recover A(l) and A' (I) 
by solving a linear algebraic system and get 



A(l) = 



C/J, 



ki - k 2 



ki 



Z(h,l) Z(k 2 ,l) 



A'{1) 



2%c[ik\k 2 



Z{k 2 ,l) Z(k u l) 



(6.4) 



ki - k 2 

On the other hand, if we only know \Z(k, l)\ without knowing its phase, then from (6.3) 
we get A(l) and |A'(/)| as 



A{1) = 



A' {If = 



cV 

h 2 — 

Ac 2 /J 2 k 2 k 2 
h 2 h 2 



k\ 



Z(hJ)\ 2 \Z(k 2 J)\ 2 
1 1 



\z(k 2 ,i)\ 2 



(6.5) 



(6.6) 



As seen from (6.2), Z(k,l) by itself contains no other information related to Q, or A. 



7. RECOVERY FROM PRESSURE AT A MICROPHONE 

Let us place a microphone at a radial distance r from the lips and measure at that 
microphone the absolute value of the pressure resulting from the glottal volume velocity 
given in (4.1). With the help of (4.2) it follows that this is equivalent to having \P(kJ + r)\ 
at hand for k G R + . From (4.17) we see that \P(k, I + r)| is an even function of k G R, 
and from (4.16) we get 



\P(k,l + r) 



Ail) 



ckfx 

A^J A(0)\F a (k)\ 2 
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k 2 + 



A' (If 



keK. 



(7.1) 



In this section we show that the data set {\P(kJ + r)| : k G R + } by itself does not 
uniquely determine any of Q(x), n(x), or A(x), and how additional data may be used for 
the unique determination. 

Theorem 7.1 The data set {\P(k,l + r)\ : k G R + , A(l), \A'(l)\} uniquely determines each 
of Q(x), rj(x), and A(x) for x G (0, /). 

PROOF: It is known (cf. (3.9) of [21]) that for any fixed a G (0, w) we have 



\F a (k)\ = k + 0(l), k^+oo. 



(7.2) 



Thus, from (7.1) we obtain 



\P(k,l + r)\= C ^^ [l + 0(l/k)], 
4nr a/ A(0) 



k — > +oo, 



or equivalent ly 



^oj=^5 lim 

47rr fc^+oo 



\P(k,l + r) 



(7.3) 



Using (7.3) in (7.1) we get 

k 2 _ 4\P(kJ + r)\ 2 



\F a (k)\ 2 4k 2 + A'(l) 2 /A(iy 



lim 



k' 2 



ken. 



(7.4) 



k^+oo \P(kJ + r)\ 2 _ 

From (4.17) and (7.4) we see that our data set uniquely determines \F a (k)\ for k G R, and 
hence, as indicated in Theorem 3.1, Q(x) is uniquely determined for x G (0, /). Next, from 
Theorem 5.1 it follows that rj(x) is also uniquely determined for x G [0,/]. Finally, from 
(7.3) we see that A(0) is also determined by our data set, and thus we recover A(x) for 
x G [0, /] uniquely by using (5.4). | 

Note that we assume that A(l) and |A'(/)| do not change with k, and hence they are 
constants. As indicated in Section 6 they can be obtained via (6.5) and (6.6) by measuring 
the absolute value of the impedance at the lips at two different frequencies. 

Theorem 7.2 The data set {\P(kJ + r)\ : k G R + , \A'(l)\/A(l)} uniquely determines each 
of Q(x) and rj(x) for x G (0,/), and it determines A(x) for x G [0,/] up to a multiplicative 
constant. 
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PROOF: From (2.4), (4.17), and (7.4) we see that \F a (k)\ for k G R is uniquely determined 
by our data set, and hence Q(x) via Theorem 3.1 and rj(x) via Theorem 5.1 are uniquely 
determined for x G (0,/). Furthermore, from (3.5) and (4.15) we see that if we multiply 
each of A(0), A{1), and |A'(/)| by the same constant, we do not change \p(l + r, t)\ or 
equivalently we do not change \P(k,l + r)\ for k G R. Thus, our data set corresponds to 
the one-parameter family for A(x), where the parameter A(0) appears as a multiplicative 
parameter in (5.4). | 

With the help of (7.4) and Theorem 7.2 we have the following conclusions. 

Corollary 7.3 Corresponding to the data set {\P(k,l + r)\ : k G R + ,A(/)}, in general 
there exists a one-parameter family for each of Q(x), rj(x), and A(x), where \A'(l)\/A(l) 
can be chosen as the parameter. 

Corollary 7.4 Corresponding to the data set {\P(k,l + r)\ : k G R + }, in general there 
exists a two-parameter family for each ofQ(x), rj(x), and A(x), where A(l) and \A'(l)\ can 
be chosen as the parameters. 

8. RECOVERY FROM THE PRESSURE AT THE LIPS 

Let us consider the recovery of A{x) for x G [0, /] from the absolute value of the 
pressure at the lips resulting from the glottal volume velocity in (4.1). From (4.2) we see 
that our data set is equivalent to {\P(k,l)\ : k G R + }. With the help of (2.4), (2.5), 
(2.7), and (4.11), we notice that \P(k, l)\ is an even function of k G R, and hence we have 
our data actually available for k G R. In this section we show that this data set uniquely 
recovers each of Q(x), rj(x), and A(x) for x G (0,/), and we outline an explicit procedure 
to determine these quantities. 

Theorem 8.1 The data set {\P(k,l) \ : k G R + } uniquely determines each ofQ(x), rj(x), 
and A(x) for x G (0, /). 
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PROOF: From (2.2), (2.4), and (4.11) we get 

CjJL \k\ 



\P(k,l)\ = 



jMt) y/A(fij\F a (k)\ 

Using (7.2) in (8.1), we obtain 

^A(0) A(l) =c/jl( lim 



ke~R. 



(8.1) 



P(M) 



and hence 



\Fa(k)\ = 



\k\ 



lim \P(k,l)\ 



k G R. 



(8.2) 



(8.3) 



\P(k,l)\ \k^+oo 

As seen from (2.4) and (8.3), by measuring the absolute value of the pressure at the lips 
for k G R + , we get \F a (k)\ for feeR. Then, by proceeding as in Section 3, we can recover 
Q(x) for x G (0,/) and the constant cot a; appearing in (4.9). Next, by proceeding as in 
Section 5, we determine r}{x) for x G (0,/). Note that y/ A(0) A(l) is uniquely determined 



from our data via (8.2). Furthermore, as seen from (5.3), we have rj(l) = y/A(l)/A(0). 
Thus, we obtain 



A(0) = 



. lim 

Tj(l) \k^+oo 



1 



P(k,l) 



and hence we get the area function uniquely via (5.4). | 

9. RECOVERY FROM THE TRANSFER FUNCTION 

The transfer function T(k, I) from the glottis to the lips is defined as 

T (k I) - V{l,t) 

and as we see from (4.1), (4.18), and (9.1), we have 



(9.1) 



T(k,l) = 



y/A(J)e- ikl 



^W)F a (-k) 
Hence, with the help of (2.4), we get 

Mi) 



|T(fc,l)l' 



A(0)\F a (k)\< 



k+ lA ' {l) 



2 , -I'C/) 2 



k G C+. 



(9.2) 



k z + 



4A(iy 



keK. 



(9.3) 
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Using (7.2) in (9.2) we obtain 



|T(fc,Z)| = ^^[l + 0(lA)], k^+oc, 



and as a result we can recover -4(0) as 

4(0) = 



4(0 



lim |T(fc,Z)| 



2 ' 



(9.4) 



Thus, from (9.2) and (9.4), with the help of (2.4), we have 



\F a (k)\ : 



lim \T(k,l)\ : 



k — >+oo 



|T(fc,Z)| : 



I A' {If 



k 2 + 



4 A(iy< 



k e R. 



Comparing (9.3) with (7.1) we see that 



|T(M)|' 



47rr 
ck[i 



\P(k,l + r) 



keK, 



and hence we have the following conclusion. 



Corollary 9.1 For each fixed r > 0, the information contained in the data set {\T(k, l)\ : 
k G R + } is equivalent to that in {\P(k,l + r)\ : k G R + }. 

In other words, measuring the absolute value of the pressure at a microphone placed 
at some distance from the lips is equivalent to measuring the absolute value of the transfer 
function from the glottis to the lips. Consequently, we have the following analogs of the 
results of Section 7. 

Corollary 9.2 The data set {|T(fc, Z)| : k G R + ,4(Z), |4'(/)|} uniquely determines each 
of Q(x), rj(x), and A(x) for x G (0, 1). 

Corollary 9.3 The data set {|T(fc, Z)| : k G R + , |4'(/)|/4(/)} uniquely determines each 
of Q(x) and n(x) for x G (0, /) and it determines A(x) for x G [0, /] up to a multiplicative 
constant. 

Corollary 9.4 Corresponding to the data set {\T(kJ)\ : k G R + ,4(/)}, in general there 
exists a one-parameter family for each ofQ(x), rj(x), and A(x), where |4'(/)|/4(/) can be 
chosen as the parameter. 
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Corollary 9.5 Corresponding to the data set {|T(/c,/)| : k G R + }, in general there exists 
a two -parameter family for each of Q(x), rj(x), and A(x), where A(l) and \A'(l)\ can be 
chosen as the parameters. 

10. RECOVERY FROM A GREEN'S FUNCTION AT THE LIPS 

In this section we show that the absolute value of a Green's function for (1.3) at the 
lips measured for k G R + enables us to uniquely construct each of Q(x), rj(x), and A(x) 
for x G (0, /). 

The Green's function at the lips can be defined [17] as the solution <&(/,£) given in 
(1.2) when the glottal volume velocity is as in (4.1). Thus, from (1.2), (4.2), and (4.11), 
we get the Green's function at the lips as 

— rku e ik ( ct - 1 ) 

G(M;f) = -7=^ • 

y/A^)F a (-k) 

Hence, with the help of (2.4) we obtain 

|G(Mit)l = y4i'£wr fceR ' (101) 

Theorem 10.1 The data set {\G(k,l;t)\ : k G R + } uniquely determines each of Q(x), 
rj(x), and A(x) for x G (0, /). 

PROOF: From (2.4) and (10.1) it follows that \G(k, I; t)\ is independent of t and is an even 
function of k on R, and hence our data can be extended from k G R + to k G R. Using 
(7.2) in (10.1) we get 

v /I W= i im |G(M; t) r (10 ' 2) 

and hence 

\k\ \G(k,l;t)\ 



\F a (k)\ lim |G(fe, 



k G R. 



Thus, we get \F a (k)\ for k G R whenever we have \G(k,l;t)\ for k G R + . Then, as in 
Section 3 we construct Q(x) for x G (0,/) and cot a. Next, as in Section 5, we construct 

21 



rj(x) for x G [0,/]. Finally, with the help of (5.4) and (10.2) we obtain 

c 2 n 2 [r](x)] 2 



A(x) 



lim |G(fc,Z;*)| 

k— >+oo 



2 • 



Thus, the proof is complete. | 



11. RECOVERY FROM THE IMPEDANCE AT THE GLOTTIS 

The impedance at the glottis is defined as 

*<*■■» (1L1) 

Using (4.1), (4.2), and (4.11) in (11.1) we get 

„ , ckuf(-k.O) 

^-4fh^y (1L2) 

From (2.4), (2.5), (2.7), and (11.2) we see that \Z(k, 0)| is an even function of k G R, and 
hence \Z(k, 0)| is known for k G R if it is known for k G R + . In this section we show that 
the information contained in {\Z(k, 0)| : k G R + } enables us to uniquely construct Q, rj, 
and A. 

Theorem 11.1 The data set {\Z(k,0)\ : k G R + } uniquely determines each of Q(x), 
rj(x), and A(x) for x G (0, /). 

PROOF: With the help of (2.4) and (11.2) we get 

17/7, nM c\k\n\f(k,0)\ 

|Z(M)I= A(0)\F a (k)\ - fcGR ' 

and hence, using (7.2) and the fact [18-20,22] that f(k,0) = 1 + 0(l/k) as k -> +oo, we 
obtain 

\Z(k,0)\ = ^[l + O(l/k)], k^+oo, 



which leads to 



A(0) = (11.4) 

1 } lim Z fc,0 1 J 
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Thus, we can write (11.3) as 



kf(k,0) 



F a (k) 



1^°)" (11.5) 



lim \Z(k,0)\ 



As seen from (11.5), the recovery from |Z(fc, 0)| for k G R + is equivalent to the recovery 
from |/(fc,0)/F a (fc)| for G R. 

Recall that we assume that the half-line Schrodinger equation with the boundary 
condition (2.1) does not have any bound states. In that case, it is known [21] that 



kf(k,0)/F a (k) is analytic in C + , continuous in C+, nonzero in C+ \ {0}, and it is ei- 
ther nonzero at k = or has a simple zero there, and 

fc/(M) 



F a (k) 



= l + 0(l/fc), £;^ooinC+. 



As a result, we can recover kf(k, 0)/F a (k) for k G C+ from its amplitude known for k G R 



fc/(M) ( i r ,io g \tf(t,o)/F a (t)\ 



exp — / eft , k G C+ 



via 

~F a (fc) ""^ V.27T 7.^"" t-k-i0+ 
Having constructed kf(k, 0)/F a (k) for k G R, we can use (3.2) and obtain A a (k) for G R. 
Next, using (3.4) we construct \F a (k)\ for k G R. Then, as in Section 3 we construct 
for x G (0, /), and as in Section 5 we construct rj(x) for x G [0, /]. Finally, A(x) is constructed 
via (5.4) after obtaining A(0) from (11.4). | 

12. RECOVERY FROM THE REFLECTANCE AT THE GLOTTIS 

The reflectance at the glottis is defined [3,4] as the ratio of the right-moving (from 
the glottis towards the lips) pressure wave to the left-moving pressure wave. From (4.11) 
that reflectance is seen to be equal to L(—k) because the Jost solution to the full-line 
Schrodinger equation has the extension 
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where we recall that L(k) is the left reflection coefficient appearing in (2.5). In this section 
we show that either of the real and imaginary parts of the reflectance at the glottis known 
for k G R + enables us to construct Q(x) uniquely for x G (0,/). On the other hand, it 
cannot uniquely determine either rj(x) or A(x). 

Theorem 12.1 Either of the data sets {Re[L(k)} : k G R+} and {Im[L(k)] : k G R+} 
uniquely determines Q(x) for x G (0,/). On the other hand, to each of these data sets, 
there correspond a one-parameter family for n(x) and a two-parameter family for A(x). 

PROOF: From (2.7) we see that Re[L(k)] is an even function of k on R and Im[L(/c)] is 
an odd function. Hence, both of our data sets can be extended from k G R + to k G R. In 
the absence of bound states, because Q(x) = for x < 0, it is known [22,24] that L(k) is 
analytic in C + , continuous in C+, and o{l/k) as k — > oo in C+. Thus, we can construct 
L(k) from either data set via the Schwarz integral formula [cf. (3.4)] as 



We can then obtain Q uniquely, for example, via [22,24] the Faddeev-Marchenko method, 
as 

Q( x ) = 2-^B r (x,0 + ), iGR, 
ax 

where B r (x,y) is obtained by solving the (right) Faddeev-Marchenko integral equation 

poo 

B T (x, y) + L(-2x + y) + / dz L(-2x + y + z) B T (x, z) = 0, iGR, y > 0. 

Jo 

It is already known [19,20,22,24] that T(k) and R(k) can be constructed from L(k) with 
the help of (2.6) and 



and hence, as we see from (5.8), rj(x) is not uniquely determined and we have the corre- 
sponding one-parameter family for rj(x) with cot a being the parameter. From (5.1) we 




k G C+. 




k G C+, 
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see that we have the corresponding two-parameter family for A(x), where the parameters 
can be chosen, for example, as A(0) and A'(0), or as A(l) and A' (I). | 

Finally, let us remark that, as indicated in (6.4), we can obtain A(l) and A' (I) if we 
know the impedance at the lips at two distinct real k values. 



13. EXAMPLES 



In this section we illustrate the theoretical results presented in the previous sections 
with some examples. 

Let us use / = 17.5 cm, c = 3.43 x 10 4 cm/sec, /j, = 1.2 x 10~ 3 gm/cm 3 , A(0) = 5 cm 2 , 
and A'(0) = -0.52 cm, and 

80(7 + 3^5) e 2VEx 



Q(x) 



(7 + 3y/E) e 2 ^ x - 2 



When Q is viewed as a potential of the full-line Schrodinger equation with support on R + , 
the corresponding scattering coefficients r(k), p(k), £(k) and the left Jost solution g\(k,x) 
introduced in Section 5 are rational functions of k, and it can be verified that 



9l (k,x) = e ikx 
r(k) = 



1 + 



k + i\fh (7 + 3>/5) e 2VEx - 2 
k(k + iV5) 2 



x > 0, 



£(k) 



(fc + i)(fc + 2i)' v ' (k + i)(k + 2iY 

(k + i)(k + 2i)(k-iV5) 
All the quantities related to (1.1), (1.3), and (1.6) can now be explicitly evaluated. For 

example, the left Jost solution g\(k, x) for x < for the full-line Schrodinger equation can 

be obtained as 



gi(k,x) 



+ 



x < 0. 



r(k) r{k) ' 

Via (4.9) we get cot a = —0.052, rj(x) can be obtained via (5.6), A(x) via (5.4), f(k, x) via 
(5.9), F a {k) via (2.3), the scattering coefficients T{k), R(k), and L(k) via (3.9) and (3.10), 
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P(k,x) via (4.11), v(x,t) via (4.18), \P(k,l + r)\ via (7.1), \P(k,l)\ via (8.1), \Z(k,l)\ via 
(6.3), \T(k,l)\ via (9.3), \G(k,l;t)\ via (10.1), |Z(fc,0)| vi a (11.3), A a (k) via (3.2). We also 
compute A(l) = 11.596 cm 2 and A' (I) = 0.681 cm. Even though all these quantities can 
be explicitly written in terms of elementary functions in closed forms, the corresponding 
expressions are too long to display here, and instead we only show some of their graphs. 




Fig. 13.1 Q(x). Fig. 13.2 A(x). Fig. 13.3 \P(k,l)\. 




Fig. 13.4 \Z(k,l)\. Fig. 13.5 |Z(fc,0)|. Fi S- 13.6 \G(k,l;t)\. 




k 



Fig. 13.7 \P(k, I + 20)\. Fig. 13.8 |T(fc,Z)|. Fi S- 13-9 L(-k). 

Let us remark that in Fig. 13.7 the absolute value of the pressure is measured at a 
distance of 20 cm from the lips. In Fig. 13.9 the real and imaginary parts of the reflectance 
are displayed, which can be told apart from Re[L(0)] = —1 and Im[L(0)] = 0. 

It is known from Section 6 that the graph in Fig. 13.4 cannot determine either of 
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Figs. 13.1 and 13.2. We know from Section 8 that the graph in Fig 13.3 uniquely determines 
the graphs of Figs. 13.1 and 13.2. As we know from Section 11, the graph in Fig 13.5 
uniquely determines the graphs of Figs. 13.1 and 13.2. From Section 10 we know that the 
graph in Fig 13.6 also uniquely determines the graphs of Figs. 13.1 and 13.2. We know 
from Section 9 that the graphs of Figs. 13.7 and 13.8 contain the same information, but 
neither is sufficient to determine uniquely either of the graphs in Figs. 13.1 and 13.2. On 
the other hand, as we have seen in Sections 7 and 9, the graphs in Figs. 13.4 and 13.7 
together, or the graphs in Figs. 13.4 and 13.8 together, uniquely determine the graphs of 
Figs. 13.1 and 13.2. We know from Section 12 that either graph in Fig. 13.9 uniquely 
constructs the graph in Fig. 13.1, but not that in Fig. 13.2. 
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